Optical model potentials involving loosely bound p-shell nuclei around 10 MeV/A 
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We present the results of a search for optical model potentials for use in the description of elastic 
scattering and transfer reactions involving stable and radioactive p-shell nuclei. This was done in 
connection with our program to use transfer reactions to obtain data for nuclear astrophysics, in 
particular for the determination of the astrophysical S17 factor for ^Be(p, 7)*B using two (^Be,*B) 
proton transfer reactions. Elastic scattering was measured using ^Li, ^"B, ^''C and ^^N projectiles on 
^Be and ^^C targets at or about E/A=10 MeV/nucleon. Woods-Saxon type optical model potentials 
were extracted and are compared with potentials obtained from a microscopic double folding model. 
Several nucleon-nucleon effective interactions were used: M3Y with zero range and finite range 
exchange term, two density dependent versions of M3Y and the effective interaction of Jeukenne, 
Lejeune and Mahaux. We find that the latter one, which has an independent imaginary part, gives 
the best description. Furthermore, we find the renormalization constant for the real part of the 
folding potential to be nearly independent of the projectile-target combination at this energy and 
that no renormalization is needed for the imaginary part. From this analysis, we are able to eliminate 
an ambiguity in optical model parameters and thus better determine the Asymptotic Normalization 
Coefficient for i°B^^B-Hp. Finally we use these results to find optical model potentials for unstable 
nuclei with emphasis on the reliability of the description they provide for peripheral proton transfer 
reactions. We discuss the uncertainty introduced by the procedure in the prediction of the DWBA 
cross sections for the (^Be,*B) reactions used in extracting the astrophysical factor Si7(0). 



I. INTRODUCTION 



Transfer reactions have been proposed as an indirect method to determine direct capture reaction rates at stellar 
temperatures for some time Recently we used the Asymptotic Normalization Coefficient (ANC) method to 

determine the cross section for the radiative proton capture process ''Be(p,7)®B at solar energies, or equivalently, the 
astrophysical factor, 817(0). The method relies on the fact that at low energies a capture reaction to a loosely bound 
state is a surface process. Its cross section is determined by the tail of the radial overlap integral between the bound 
state wave function of the final nucleus and those of the initial colliding nuclei. This overlap integral is asymptotically 
proportional to a well known Whittaker function, and therefore the knowledge of its asymptotic normalization alone 
determines the cross section. This asymptotic normalization, in turn, can be determined from the measurement of a 
transfer reaction involving the same vertex, provided that this second reaction is also peripheral. In particular, we 
determined 817(0) from measurements of the ANC for the *B— >^Be -|- p system utilizing the proton transfer reactions 
^'^B(^Be,^B)^Be |Q] and ^''N(^Be, '*B)^'^C 1^, at energies where the proton transfer process is peripheral. Determining 
the ANCs from transfer reactions involves distorted wave Born approximation (DWBA) calculations, and therefore 
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good, reliable optical potentials are needed. In particular, good optical model potentials are needed in both the initial 
and the final channels involving ^Be and ®B radioactive nuclei in each reaction in order to compute the DWBA proton 
transfer cross sections. Because elastic scattering data with these projectiles cannot easily be obtained and/or are not 
precise enough to extract reliable and unambiguous optical potentials, we have studied elastic scattering for several 
combinations of p-shell nuclei at energies close to those appearing in the reactions of interest. We then determine 
procedures to extract optical model potentials for the cases involving radioactive partners. 

Angular distributions up to the nuclear rainbow angle were measured in seven experiments using ^Li, ^°B, ^'^C and 
^'*N projectiles on ^Be and ^'^C targets at bombarding energies at or around E/A=10 MeV/nucleon. They were fitted 
with phenomenological potentials with volume Woods-Saxon real and imaginary terms. The phenomenological optical 
model potentials found for all systems were then compared with the potentials calculated with microscopic double 
folding procedures, using six effective nucleon-nucleon interactions: M3Y with zero range and finite exchange term, 
the density dependent M3Y interaction in the forms extracted recently by Khoa et al. ||] (BDM3Y1, BDM3Y3) and 
the interaction of Jeukenne, Lejeune and Mahaux (JLM) |^ in two versions. For the calculations, one-body densities 
were obtained in a standard spherical Hartree-Fock calculation using the density functional of Beiner and Lombard 
with a slight modification of the surface term in order to fit the experimental binding energy for each nucleus. 
These densities were used in the double folding procedure. The renormalization coefhcients needed for the analysis 
of elastic data with these double folding potentials are extracted and discussed here. 

In two cases the elastic scattering data were also used to extract the parameters of the optical model potentials that 
were needed for DWBA calculations to determine the ANCs for the ^''B— >^Be -I- p and ^^N ^^■^C-l- p systems from 
the ^Be(i°B,9Be)i°B @ and ^3c(14n,i3c)14n reactions, respectively, and were included in those publications. 
These measurements also allowed us to test the results of our folding procedures for the proton transfer reactions 
by comparing the results of the DWBA calculations that we obtain from the double folding model potentials with 
those we obtain using the phenomenological potentials. From the analysis presented below, we are able to eliminate 
the uncertainty in the optical potentials found in Ref. This results in a better determination of the ANC for 
iog^9gg_l_p -y^g chose the elastic scattering of '^Li on both targets as being close to what we expect for the scattering 
of ^Be. Both ^Li and ^Be nuclei have low binding energies and large break-up probabilities. The elastic scattering of 
^^C on the ^Be target was studied as this is close to the exit channel ^'^C-|-*B of the second proton transfer reaction 
"N(7Be,8B)i3C. 

The double folding procedure and the parameters extracted for the effective nucleon-nucleon interaction from the 
present experimental data were also checked for other projectile-target combinations in the same mass and energy 
region for which data exist in the literature. The results were in agreement with those found from the study of the seven 
cases described here. Using the techniques developed here, we calculated the optical model parameters required to 
extract the ANC for the ^Be-|-p^^B system, and consequently 817(0), from the studies of the ^*'B(^Be,^B)^Be |^] and 
^'*N(^Be,*B)^'^C H] reactions. Section || presents the experiments and the procedures used in the data reduction. We 
extract the Woods-Saxon optical model potentials from fits to the elastic scattering data in Section |l|, and compare 
them with those microscopically calculated in Section [V. In sections ^ and VI, we describe a global optical potential 
for interactions involving l oosel y bound p-shell nuclei and its implications for the ^°B(^Be,*B)^Be and ^'*N(^Be,*B)^^C 
reactions. Finally, Section VII presents the conclusions. 



II. THE EXPERIMENTS 



The experiments were carried out using beams from the Texas A&M University K500 superconducting cyclotron. 
A list of the measurements is given in Table ^. The experimental setup and the data reduction procedures were 
similar to those already described in Ref. The Multipole-Dipole-Multipole (MDM) magnetic spectrometer |Tl| ] 
was used to analyze the scattered particles and the reaction products. The beams were prepared using the beam 
analysis system [n2j, which allows for the control of the energy and angular spread of the beam. Self-supported ^Be 
and ^^C targets, about 200-300 /ig/cm^ thick, obtained by evaporation, were placed perpendicular to the beam in the 
sliding-seal target chamber of the MDM. The magnetic field of the spectrometer was set to transport fully stripped 
ions to its focal plane, where they were observed in the modified Oxford detector There, the position of the 

particles along the dispersive direction was measured with resistive wires at four different depths within the detector, 
separated by about 16 cm each. For particle identification we used the specific energy loss measured in the ionization 
chamber and the residual energy measured in a NE102A plastic scintillator located in air, just behind the output 
window of the detector. The input and output windows of the detector were made of 1.8 and 7.2 mg/cm^ thick 
Kapton foils, respectively. The ionization chamber was filled with pure isobutane at pressures between 30 — 50 torr. 
The entire horizontal acceptance of the spectrometer, = ±2°, and a restricted vertical opening, A(f) = ±0.5°, were 
used in these measurements. Raytracing was used to reconstruct the scattering angle in the analysis of the data. For 
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this purpose, in addition to RAYTRACE jlj] calculations, angle calibration data were obtained at several angles by 
using an angle mask consisting of five openings of 69=0.1°, located at -1.6°, -0.8°, 0°, -1-0.8° and -1-1.6° relative to the 
central angle of the spectrometer. Moving the spectrometer from Oiab = 4° to 54° we covered the angular ranges listed 
in Table |. Typically we moved the spectrometer by 2° or 3° at a time, allowing for an angle overlap that provided 
a self-consistency check of the data at all angles. Normalization of the data was done using current integration in a 
Faraday cup. Focal plane reconstruction was done at each angle using the position measured with the signals in the 
wire nearest to the focal plane and using the detector angle obtained from the position measured at two of the four 
wires. The angular range of A6 = 4° covered by the acceptance slit was divided into 8 bins, resulting in 8 points in 
the angular distribution being measured simultaneously, with each integrating over 69iab = 0.5°. 

The measurements with the angle mask showed that the resolution in the scattering angle (laboratory) was AOres — 
0.20° — 0.25° full-width at half maximum (FWHM). This includes a contribution from the angular spread of the beam 
of about 0.1°. The best energy resolution obtained at forward angles was 150 keV FWHM for on ^Be, 230 keV 
for i^'N on and 150 and 220 keV for the scattering of "^Li at 63 MeV and 130 MeV, respectively. It degraded 
as we advanced to larger angles due to the large kinematic factor, k = coupled with the finite angular spread 
in the beam. However, it never degraded our ability to isolate the elastic peak, even in the most demanding case of 
the ''Li experiments where the first excited state of the projectile is only 477 keV away. A sample spectrum taken 
in one of these most demanding cases, ^Li (130 MeV) -I- ®Be, is presented in Figure 1 where the good separation is 
clear. The active length of the focal plane allowed us to cover a total excitation energy of about 7 MeV, centered 
around the elastic peak. Thus we were able to measure inelastic scattering to the lowest excited states of the 
projectile-target systems at the same time. These inelastic scattering data were used as additional information to 
check the experimental procedures, and in a few cases we compared the inelastic transition strength obtained in these 
experiments (deformation parameters or deformation lengths) with those available in literature. 

During the experiments, particular emphasis was placed on obtaining accurate absolute values for the cross sections, 
and therefore target thickness and charge collection factors were determined by a two-target method and by normaliza- 
tion to Rutherford scattering at forward-most angles, as described in Combining the results of these independent 
determinations, we conclude that we have an overall normalization accuracy of 7% for the absolute values of the cross 
sections for both the elastic scattering ^Be(^*'B,^"B)^Be data and for the proton transfer ^Be(^"B,^Be)^°B data j|] 
and for the elastic scattering 13c(14n,i4n)13c and proton transfer reaction i3c(i''N,i3c)i4]y rj.^^ normalization 
for the absolute values of the cross section was made in the ^^C (130 MeV) -|- ^Be case only using the nine most 
forward angle data points in the angular distribution and is estimated to be accurate to 10%. For the experiments 
with the ^Li beam, we also determined the target thickness by measuring the energy loss of alpha particles from a 
^^^Th source and the accuracy in normalization is 9%. 



III. PHENOMENOLOGICAL ANALYSIS 



The elastic scattering data have been fitted using the code OPTIMINIX in a standard optical model analysis 
using volume Woods-Saxon form factors with the standard notation: 



U(r) = -(Vfv(r)+iWfw(r)) 



(1) 



where 



fx(r) = 



1 + exp 



(2) 



and x=V,W stands for the real and imaginary parts of the potentials, respectively. Fits using the codes PTOLEMY 
Ip^ and ECIS [0 gave similar results. Only the central components have been included in the optical potential, 
since vector and higher rank tensor spin-orbit couplings have little or no influence on the cross sections. 

In Figure || we present the angular distributions measured for ^"B+^Be at E(^°B)=100 MeV, ^''N-I-^^C at 
E(i'*N) = 162 MeV and ^^C+^Be at E(i3C)=130 MeV and in Figure | those for ^U+^Be,^^C at E(^Li)=63 and 
130 MeV. All angular distributions display typical patterns for elastic scattering, dominated by strong absorption 
with Fraunhoffer oscillations of large amplitude around the crossing point, followed by less developed structures at 
larger angles. The Fraunhoffer oscillations are expected for such systems, due to the small Sommerfeld parameter 
77 ~ 1. The curves are the fits to the data. Inspecting Figures g and ^, one observes that all potentials predict 
somewhat deeper minima than the data show. This effect is partially attributable to the finite angular resolution. 
The rest is probably due to the incoherent contribution of the substantial quadrupole moment of some of the partners 
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(like ^Be). During the normalization procedure, the theoretical curves were convoluted with the experimental angular 
resolution and binning, using the code ECIS, but this was found to have no influence on the fits. 

The optical model parameters extracted for all seven cases studied here are presented in Table ||. In addition to 
the depth, reduced radius and diffuseness for the real and imaginary parts of the potential, the table gives the chi- 
square for the fit (x^), the total reaction cross section (cfl), the values of the volume integrals per pair of interacting 
nucleons for the real (J\/) and imaginary parts {3w) of the potential, respectively, and the root- mean-square radii of 
the real (Ry) and imaginary (Rh/) potentials. The parameters were obtained by griding the initial strength of the 
real potential in small steps in the range V=50-250 MeV in order to identify the local minima and then searching 
for minima on all parameters with no constraints. In this way 3-4 families of potentials have been found for each 
case. Usually, a characteristic jump of 50-70 MeV fm'^ in the volume integral of the real part of the optical potential 
serves to identify these potentials as discrete members of a sequence of potentials which give a comparable description 
of the data. The members of each family of potentials are connected by the well known continuous Igo ambiguity: 
Vexp{Rv/av) — const. 18|. This arises since, due to the strong absorption, the cross section is sensitive only to 
the tail of the potential. The Igo potentials of the same family have practically the same volume integral and the 
same radius and therefore the discrete families of potentials can be identified by the values of these parameters. The 
absorption is seen to be independent of the strength and shape of the real part of the optical potential and, as a 
consequence, the reaction cross section is more or less constant along the sequence in each case. Also we notice that 
generally the radii of the imaginary potentials are about 20% larger than those for the real potentials, in agreement 
with previous observations |19|| . 

For the ^°B-|-^Be case, it appears that potential 1, which has the smallest chi-square, provides the most realistic 
description, and potential 3 can be rejected. In the angular range covered, the prediction of potential 2 for the elastic 
scattering differs from that of potential 1 primarily in the depths of its minima. We reached the same conclusion 
from the comparison of the DWBA calculations for the proton transfer reaction ^Be(^"B,^Be)^°B studied in the same 
experiment: whereas potentials 1 and 2 give a very good description of the shape of the angular distribution and 
similar absolute values, potential 3 predicts a reaction cross section which is far too small To further clarify the 
features of the angular distribution we have performed a near-far decomposition of the scattering amplitude, with one 
potential taken as the reference potential. Shown in Figure ^ are the cross-sections due to the near-side and far-side 
components of the total scattering amplitude. Around the crossing angle of 9c.m. = 16°, the strong interference 
between the near- and far-side components results in Fraunhoffer oscillations of large amplitude. Beyond this angle, 
the near-side component drops rapidly and the angular distribution is dominated by the far-side component, which 
falls off smoothly. No significant structure could be identified in this region. This structureless behavior can be 
understood qualitatively in terms of the transparency for the low partial waves implied by the refractive potential. 
In the case of "'^"B -I- ^Be, the rainbow partial wave has Ir — 16 and the associated scattering matrix elements are of 
the order \Si \ ~ 3 • 10"'^ for I << Ir. Thus, the refractive nature of the potential is sufficient to allow the interference 
between waves with I < 16 and higher ones to produce the smooth behavior. Comparison of the potential elastic 
scattering branch (forward angles) and of the elastic proton exchange branch (backward angles) in Figure 3 of Ref. 
1^ shows that the interference between these two mechanisms has no sizable effect in the angle range considered here 
and was not considered in the analysis. The data for the ^■^C (130 MeV) -|- ^Be experiment show similar features, 
with 3 families of potentials found. 

For the -I- ^^C system the volume integral and radius for the absorptive part of the optical potential seems to be 
independent of the real potential, resulting in a constant total reaction cross section along the sequence with an average 
value of 1535 mb. The optical model total reaction cross section is consistent with the experimental value measured by 
DiGregorio et al. at 161.3 MeV, a = 1463 ±100 mb |20|. All potentials give reasonable x^, but potential 1 listed in the 
table gives the smallest value and is the only one that fits the data at the largest angles. This potential has a volume 
integral per pair of interacting nucleons close to that which we found for ^°B-|-^Be elastic scattering at similar velocities. 
Potential 1 was adopted for the DWBA calculation of the proton exchange process ^'^C(^^N,^'^C)^'^N as described in 
Ref. jl^, while the others were used to estimate the uncertainty due to the choice of optical model parameters. Similar 
insight on the relative role played by the refractive and absorptive parts of the optical potential may be obtained 
from the far side- near side decomposition of the scattering amplitude corresponding to potential 1 which is presented 
in Figure ^p. The far-side component is represented by the dashed line and the near-side component by the dotted 
line and their coherent sum by the continuous line. For angles around the crossing where the two components have 
comparable amplitude and strongly interfere, a typical Fraunhofer diffraction pattern emerges with large amplitude 
oscillations equally spaced by A6' = Tr/lg = 7r/30, where Ig is the grazing angular momentum. Beyond this angle, the 
near side component is completely damped by the strong absorption and we are left with the far side exponential tail 
that is characteristic of far-side dominance. No significant structure could be identified up to the nuclear rainbow 
angle, which in this case is Or — 83°. Similar to the case of the ^°B-|-^Be experiment, our measurements show that 
we do not have interference effects between the potential scattering predominant at forward angles and the elastic 
proton exchange predominant at backward angles (see Figure 2 of Ref. [pO|) in the angular range considered. 
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The potentials found in the phenomenological WS analysis of ''Li scattering are given in Table ||. A similar result 
emerges, with discrete ambiguities represented by up to 4 families found in each case. Similar values are found for 
the volume integrals for the real and imaginary parts as for the rest of the cases studied above. We notice however 
that the reduced radii ry^w are small and the diffusivities ayw of the potentials are unusually large. This agrees 
with findings in other analyses for such light systems. Figure ^ shows the far-side, near-side decomposition for the 
''LiH-^'^C system at 63 MeV, with conclusions similar with those for the cases discussed above. 

It is interesting to note that, for all but one of the cases shown in Table ||, the first of the potentials always has a 
similar volume integral for the real part: Jy ~ 220 MeV fm^, and that the imaginary potentials are independent of 
the real part, predicting the same total cross sections. 

As mentioned above, the spin dependent components of the optical potential have been omitted. In the absence of 
any polarization data, exploratory calculations for the ^'^B + ^Be system, using the same vector spin-orbit term as 
for ^Li + ^^C at E=156 MeV, did not result in any noticeable effects on the elastic cross section in the measured 
angular range. For several of the cases studied here, we also did a Fourier-Bessel analysis of the data, similar to that 
in Ref. 123], and found that the phenomenological Woods-Saxon shapes assumed in Eq. 2 are adequate. 



IV. FOLDING MODEL ANALYSIS 



In addition to the analysis with Woods-Saxon type potentials, the data have been analyzed in the framework of a 
semi-microscopic folding model. We followed a Hartree-Fock procedure to obtain the densities in the two partners, 
then used double folding with known nucleon-nucleon interactions. The wave functions and the densities for all 
nuclei involved were obtained in a standard spherical Hartree-Fock calculation using the energy density functional of 
Beiner and Lombard j|] . This functional describes nuclear matter and the bulk properties of finite nuclei well. In the 
calculations, the parameters of the surface terms were adjusted slightly in order to reproduce the experimental total 
binding energy. This is an important constraint on the calculation, especially for nuclei with small separation energies 
such as ^Be and ''Li. Usually this correction amounts to a few percent with respect to the original parameters and 
substantially improves the description of the single particle levels close to the Fermi level. The calculated binding 



energies and the rms radii that were obtained are given in Table [II and compared with the experimental ones. A 
similar procedure has been used by Hoshino et al. to describe the structure of the ^^Be nucleus. 

In the double folding procedure, we used a number of G-matrix effective nucleon-nucleon (NN) interactions. The 
first one is the well known M3Y interaction. Recall that the nucleus-nucleus potential in the double folding model is 
given by 



VfoidiR) = jdridr2pi{ri)p2{r2)veff{ri+R-r2) (3) 

where pi, 2 are the single particle densities, and the interaction operator is of the form 

«e//(r-) = voir) + Pi%Vex{r) 

where the direct and exchange parts are averaged over spin-isospin states and Pf 2 is the knock-on exchange operator 
in coordinate space. We assumed, as usual, that the one nucleon exchange knock-on term, which involves the exchange 
between the interacting nucleons, is dominant with respect to all other exchange contributions. The parameters for 
the direct and exchange components of M3Y were taken from Ref. |2^ . In the standard version , the isoscalar 
component of the interaction consists of a finite range direct term, supplemented by an energy dependent zero range 
pseudo-potential which simulates well the one nucleon knock-on contribution to the interaction. The small isovector 
component of the interaction has also been included in the calculation and the corresponding results are denoted by 
M3Y/ZR throughout the paper. A finite range version of the M3Y interaction was also used for some of the systems 
analyzed in this paper. The lack of any explicit density dependence in the effective interaction results in potentials 
that are too deep in the interior to reproduce correctly the rainbow features at large angles observed e.g. in alpha 
scattering at higher energy ||2^ . This can be corrected by making the effective NN interaction depend upon the density 
of the nuclear matter in which the interacting nucleons are immersed. The requirement that nuclear matter saturate 
ensures that this density dependence reduces the strength of the interaction as the density increases, weakening the 
folding potential in the interior while leaving the surface values practically unchanged. For our purpose we adopted 
more recent interactions called BDM3Yn(n=l and 3) which have been shown to give a good description of light ion 
scattering in a wide range of incident energies Q . These interactions are based upon a G matrix derived from the 
Reid soft-core NN potential. They incorporate a linear (n=l) or cubic (n=3) density dependence with parameters 
adjusted to give saturation in nuclear matter at the correct density and binding energy. The two interactions give 
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very different compressibilities for nuclear matter {K = 230 MeV for n = 1 and K = 475 MeV for n = 3) covering a 
broad range of equations of state. We note that at present Koa = 231 ± 5 MeV has experimental support ||2^, which 
would indicate a preference for BDM3Y1 in the description of heavy ion elastic scattering. 

Also, we have used the G- matrix interaction of Jeukenne, Lejeune and Mahaux (JLM) [Q, which is obtained in a 
Brueckner-Hartree-Fock (BHF) approximation from the Reid soft-core nucleon-nucleon potential. This interaction is 
complex, energy and density-dependent and, therefore, provides simultaneously both real and imaginary parts of the 
optical potential. The interaction has been considered recently by Bauge, Delaroche and Girod pO| in an extensive 
study of nucleon scattering on a wide range of target masses and incident energies. Some shortcomings of the original 
interaction were also corrected in Ref. |^^. For completeness, we describe below the main steps in the derivation of 
our potentials, taking into account the improvements recommended in Ref. |30| . 

The optical potential for a nucleon of energy E traversing nuclear matter of density p is written as: 

Unm{p, E) ^ Vo{p, E) + arVii^p, E) + i[Wo{p. E) + aTWi{p, E)] (4) 

where a = {pn — Pp)/{pn + Pp) and r — ±1 for neutrons and protons, respectively. Explicit expressions for various 
terms are : 

3 

Vo{p,E)=Y,a,,p'E^-^ (5) 



Vi{p,E)^-m{p,E) (7) 



TO 



TO 



Wi{p,E)^ -^N{p,E) (8) 



TO 



The matrix coefficients a^ , dij, the Fermi energy e_F and the BHF expression of the auxiliary function N{p,E) are 
given in pO{ . The quantities rh/m and fh/m are the k mass and the E mass, respectively, and represent a measure 
of the true nonlocality and the true energy dependence of the optical potential. 
In applications for heavy ions we interpret the quantity 

vo{p, E) = (Fo(p, E) + zWo(p, E))/p (9) 

as the (complex) isoscalar, density- and energy-dependent NN effective interaction. The heavy ion potential is given 
then by the folding integral: 

V{R)^ J drldr2Pi{ri)p2{r2)vo{p,E)S{s) (10) 

with s = 7'1 + i? — r2 . Similarly, the quantity 

vi{p,E) - iViip,E)+tWiip,E))/p (11) 

is interpreted as the (complex) isovector, density- and energy-dependent NN effective interaction. The corresponding 
heavy ion potentials are obtained from a folding integral similar to that in Eq. (10), replacing vq by vi and the single 
particle densities pi^2 by the isovector densities (/5„ — Pp)ifi. Usually such terms have little influence in the total 
optical potential because the isovector densities are small for normal nuclei in the p-shell; however we have included 
them in the analysis since such terms can have some importance in the case of loosely bound nuclei with very different 
proton and neutron single particle densities. Two approximations for the local density have been used. The first of 
them reads: 

/ - - \ 1/2 

P= Pi(^"! + J)P2(r2-^) (12) 
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which amounts to an estimate of the local density as the geometric average of the individual single particle densities, 
each of them evaluated at the mid distance between the interacting nucleons. This approximation has been used 
by Campi and Sprung in Hartree-Fock calculations with density-dependent forces. With this approximation, 
the local density never exceeds the saturation value for nuclear matter density po- We remind the reader that the 
JLM effective interaction is defined only for density values satisfying p < pq. Potentials obtained with the above 
approximation are labeled below as JLM(l). The second approximation for the local density uses the arithmetic 
average of the individual densities: 

^ Pi{rl + ^) + P2{r2 - '^)] (13) 



A similar approximation was used by the authors of Ref. in their derivation of the density-dependent version of 
M3Y, except for the factor of one-half in front of the parentheses which is introduced here in order to be consistent 
with the assumptions of the JLM model. Potentials calculated with this approximation are denoted by JLM(2). It 
has been shown by the authors of Refs. and |3^] that the local density approximation is substantially improved by 
replacing the 6 function in integrals of the type (10) by finite range form factors of the gaussian shape: 

,i^-{j^)'e-'"' (14) 

Since the finite range form factors are normalized to one, the volume integrals of the folding potentials are not affected, 
only the rms radii are increased, depending on the values one chooses for the range parameter t. Our phenomenological 
analysis shows clearly that the bulk of the elastic scattering experimental data require larger radii for the absorptive 
part of the heavy ion optical potentials as compared to the real part. Extensive numerical calculations with both 
versions of the JLM interaction showed that optimum values for the range parameters are Ir = 1.2 fm and tj = 1.75 fm. 
A similar need for different radii of the imaginary and real parts of the optical potential has been emphasized recently 
by Satchler and Khoa . Of course slightly improved fits could be obtained in each individual case by varying also 
the range parameters around these values. For example, our ^Li data were better fitted with a larger tj. However, 
finding such variations for individual data sets goes beyond the purpose of the present paper. 

It is known that p-shell nuclei elastic scattering, some of which involve loosely bound nuclei, cannot be described 
successfully without a substantial renormalization of the folding form factor . The strong coupling with breakup 
and neutron transfer channels is responsible for such an effect. The usual procedure to simulate the repulsive effect 
of the real part of the dynamic polarization potential |Q arising from such coupling is to introduce a multiplicative 
constant for the real folding form factor. In the folding model with real effective interactions the absorption is 
accounted for phenomenologically by adding an imaginary potential of the same shape as the real part: 

U(r) = (Nv -MNw)Vfo;d(r-) (15) 

whereas for the cases when the effective interaction also has an imaginary component, the renormalization is: 

U(r) - NvVfoid{r)+iNyjWfoid{r). (16) 

The resulting potentials differ from the Woods-Saxon shape at small distances, but can be easily fitted with such 
forms in their surface region. We reanalyzed all our elastic scattering data using double folding potentials obtained 
with the six effective interactions outlined above. The renormalization constants Ny and Nw were further adjusted 
to fit the elastic scattering data using Eq. |l5| in the case of M3Y and BDM3Y forces and Eq. ^ for the two versions 
of the JLM interactions. The results of the fits are shown in Figures ^ and and the parameters are displayed in 



Table [V. In general, fits of reasonable quality were obtained with all interactions. However, the JLM interaction 
not only gives the best fits as compared to the other interactions, but also provides renormalization constants with 
minimal dispersion for all projectile-target combinations considered. This indicates that the mass dependence of the 
optical potential is properly taken into account by these effective interactions through the density dependence. As a 
rule, all folding potentials need a substantial renormalization for the real part of the optical potential, emphasizing 
that the dynamic polarization potential plays an important role for p-shell nuclei elastic scattering at low energies. 
Density-dependent effects, such as those taken into account by BDM3Y forces, lead only to a slight increase in the 
real normalization constant Ny as compared to the original density-independent interaction M3Y/ZR, suggesting a 
need for a stronger density dependence at the potential surface. Inspecting Figures ^ and ^ one sees that is hard to 
distinguish between the two versions of this force since both of them give a comparable description of the data. This 
is likely a consequence of the fact that the present data give information on the optical potential in a limited spatial 
region centered around the strong absorption radius where the two do not differ much. 
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V. EXTRACTING A GLOBAL OPTICAL POTENTIAL 



The analysis done as described in the previous section leads us to the conclusion that we can find a way to 
predict optical model potentials with some reliability. As already noted before, the situation is complicated by the 
fact that the nuclei involved are loosely bound and we expect to have important effects from break-up channels. 
Satchler and Love ||l^ concluded earlier that the renormalization of the real part of the double folding potentials is 
considerable, particularly for loosely bound nuclei where break-up is important. The energies studied here, around 
10 MeV/nucleon, are known to lead to sizable effects due to the dynamic polarization contribution to the optical 
potential This is most likely the explanation behind the need for a substantially reduced real well depth. The 
renormalization coefficients are presented in Table ^ and in Figure ^ for both the real and imaginary part of the 
potentials. If we compare the results for the same nucleon-nucleon interaction, we see that similar renormalization 
constants are obtained for all systems when at least one of the participating nuclei is weakly bound. In particular 
when density-dependent effective interactions (JLM, BDM3Y1, BDM3Y3) are used, the renormalization constants 
are very stable, with a standard deviation of a few percent around the average value. This suggests that one can 
indeed obtain the optical model potentials for pairs of projectile-target nuclei for which data are not available, or 
are scarce, by using a folding procedure to obtain the geometrical parameters and the renormalization constants 
extracted above. Studies like the one comparing the scattering of ^Li and ^^Li lead to a similar conclusion, and 
show that the energy dependence in the potential is smooth and rather weak p6[ |. Furthermore, the renormalization 
factors that we find here are comparable to those found for ^Li-f^^C near this energy when the M3Y and JLM 
interactions are used . In a few cases the finite range version of the exchange term in the M3Y interaction was also 
checked but the results were not improved over those obtained with the zero range version of it. Given our suspicion 
that the localization procedure used to obtain these finite range calculations might not work properly in very light 



nuclei, we do not discuss the results here, but they are included in Table IV. From all six effective nucleon-nucleon 



interactions used above, we favor the one denoted JLM(l) because it gives a slightly better fit than the others and 
the renormalization coefficients have the smallest spread around the average value (last rows in Table IV). In contrast 
to the real potential, no renormalization is needed for the imaginary part of the calculated potential, a sign that 
the imaginary part of the effective interaction and its density dependence are well accounted for. There might be a 
remaining slight dependence of the renormalization on energy, as found in other studies, but our data are insufficient to 
extract a definite conclusion on this dependence. However, it seems likely that most of the energy dependence is taken 
care of by the energy dependence of the effective interaction and by the density dependence used in the calculations. 
We also checked our double folding procedure on other systems than those mentioned above, and included the results 



in Table IV. Whereas we obtain very good fits to the data over a large mass and energy range, thus confirming 



the appropriateness of the JLM(l) effective interaction and of the smear function and ranges used in Eq. 14, the 
resulting renormalization coefficients, when using alpha particles for example, differ from those for the p-shell nuclei 
studied here and point to the conclusion that the present coefficients have only a local applicability. Analysis of alpha 
scattering of up to 60 MeV/nucleon on stable targets, lead to renormalization coefficients for the real part about 
a factor two larger. This is surely a reflection of the differences between the very well bound ^He nuclei and the 
loosely bound partners studied here. In order to obtain more complete information on the renormalization constants, 
we have included in our analysis two angular distributions involving the elastic scattering of another loosely bound 
p-shell nucleus ^Li on light targets at 16 MeV/A and ^Li-|-^^C at two energies |3^. The volume integrals of the 
renormalized double folded potentials agree with the volume integrals of the first of the phenomenological potentials 
found and suggest that the phenomenological potentials with Jy « 220 MeV fm^ give the most realistic description. 

Data in Figure ^ and Table |^ show that the renormalization coefficient of the real potential calculated with the 
JLM(l) interaction is somewhat higher for ^^N-|--'^'^C than the average of the remaining 6 cases measured here. This is 
the only projectile-target combination where both nuclei are well bound, thus we should expect a smaller contribution 
from the polarization potential. The averages and standard deviations for all 7 cases are Nv=0.378±0.034 (or 9%) 
and Nw=1.004±0.135 (13%), respectively. Excluding the i-^NH-^^c system we find the averages Nv=0.366±0.014 (or 
4%) and Nw=1.000±0.087 (9%). We see that the value of the renormalization coefficient is very stable. We suspect 
that a large part of the spread around the average of the renormalization coefficient for the imaginary potential comes 
from the uncertainties in the absolute normalization of our data. The real part of the potential, which is fixed mostly 
by the position of the oscillations in the angular distributions, is less sensitive to this absolute normalization. This is 
the reason we exclude the data of other groups (lower part of Table IV) from the present averaging procedure (last 
row in Table |v|) . We note that, deformation, which is important in some p-shell nuclei, is not included in any way 
in our calculations, due to the use of spherical Hartree-Fock density distributions. 

Further, we checked to see to what extent the double folding potentials, renormalized to fit the elastic scattering 
data, give the same results as the phenomenological potentials when used to calculate the cross sections for the 
proton transfer reactions. For the reaction ^Be(^"B,^Be)^"B we found that the cross section calculated with the 
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JLM(l) potential, renormalized as above, differs by less than 1% (integral over the angles from 0° to 45° ) from 
that calculated using the phenomenological potential 1. Furthermore, the double folding potential used has the same 
volume integrals as the phenomenological potential 1. In Ref. ^ we left open the choice of the potential we use 
to extract the value of the ANC for the system ^Be+p^^°B, and two slightly different values were extracted using 
potentials 1 and 2. The present study and the calculations made with the double folding potential indicate that 
we can select potential 1 as the only potential, and that the value C^=4. 91(19) fm~^ is a better choice than the 
weighted average C^=5. 06(46) fm^^ given previously. For the reaction ^^C(^^N,^^C)^''N we found that the value of 
the cross section calculated using the double folding potential JLM(1) varies at any angle between 9cm = 0° — 35° 
by less than 2% from that calculated using potential 1 in Table 0, and its integral over the same angular range 
does not vary at all. This is easy to understand given the fact that the surface part of the nuclear potential is the 
contributing factor in the description of both the elastic scattering and the transfer reaction. Previously we found 
that the calculated cross sections for the proton transfer reaction ^'^C(^'*N,^'^C)^''N differ by about 2% between any 
consecutive phenomenological potential families in Table ||, as described in Ref. |jl^ . The present verification increases 
our confidence in using the double folding procedure for the description of the transfer reactions. 



VI. OPTICAL POTENTIALS FOR i°B(^BE,^B)'^BE AND ^*N('^BE,*B)i3C REACTIONS 

Using the procedure outlined above, the JLM(l) effective nucleon-nucleon interaction and the average renormal- 
ization coefficients extracted, we calculated the optical model potentials needed in the analysis of the proton transfer 
reactions i''B('^Be,^B)^Be and "N('^Be,^B)i3C at E('^Be)=84 MeV, which were the original motivation for the present 
study. The systems involve the radioactive ^Be and ^B nuclei, both loosely bound and with important clusterization 
in their ground states. This made us treat the Hartree-Fock densities carefully, forcing the calculations to reproduce 
the correct binding energies through slight modifications in the surface term as stated before. Furthermore, in the 
final calculations we imposed self-consistency by requiring that the tail of the density distribution of *B have the 
asymptotic behavior given by the ANC extracted from the transfer reaction (but not predicted by HF). This pro- 
duced changes in the potentials only at distances larger than 8 fm, as shown in Figure ^, and did not introduce any 
substantial change in the calculated proton transfer cross sections. The optical model potentials obtained reproduced 
very well the measured angular distributions for the elastic scattering of ^Be on the ^°B (Ref. Q) and ^''N (Ref. 
[^) targets without any need for further adjustments (in both cases elastic scattering was actually calculated using 
JLM(l) potentials not only on the main component of the target, ^"B and ^^N, respectively, but also on the ^^O and 
^■^C nuclei present as impurities in the "'^^B and in the melamine target, respectively). This provided confidence for 
using the extracted potentials for the description of the angular distributions of the proton transfer reactions. Again, 
the shape of the measured angular distributions are very well reproduced, as seen for "N(^Be,8B)i3c in Figure |. 
In turn, the calculated cross sections were used to extract the Asymptotic Normalization Coefficient for the system 
^B^'^Be-l-p and, consequently, the astrophysical factor 817(0) reported in Refs. |^,||. The potentials obtained for 
these nucleus-nucleus combinations involving radioactive ^Be and *B are not exactly of Woods-Saxon shape, but can 
be approximated in the region of their surface by Woods-Saxon potentials. In Table ^ we give the parameters found 
by fitting the range of radii r=2-12 fm. 

In order to estimate the uncertainty in the ANCs due to the optical model potentials, we consider that the standard 
deviations of the normalization coefficients (5Ny=0.014 and 6Nw— 0-087 give a good measure of the uncertainty with 
which we can find the depths of the real and imaginary potentials wells, respectively. By the choice of the systems 
considered here, we span a good range of p-shell nuclei, averaging properties similar to those of radioactive ones in 
terms of mass, separation energy, structure of the ground states, incident energies, number of open reaction channels, 
etc. We used these standard deviations around the average value of the renormalization coefficients to evaluate the 
uncertainty in extracting the ANCs. The uncertainties arise through the DWBA calculations of the transfer reaction 
cross section. We took the geometry as given by the double folding procedure and determined the variation of the 
calculated proton transfer cross section integrated over the angular range relevant in the experiments. The potential 
depths were varied from Ny-^Ny to Ny-|-(5Ny for the real part and from Nvi/-JNw to Nvv-|-(5Nvi/ for the imaginary part 
for the entrance and exit channels independently and the resulting variations were added in quadrature to estimate 
the relative uncertainty in the DWBA calculations. With this procedure we found a 7.5% uncertainty in the calculated 
^°BCBe,^BfBc transfer cross section due to DWBA calculations. The same procedure gave an estimate of 7.7% for 
the uncertainty due to DWBA calculations for the ^'*N(^Be,*B)^'^C reaction. Most of the contribution comes from the 
uncertainty in the imaginary renormalization coefficient (7.5%), while the real one contributes only about a quarter of 
that (2%). Note that in varying separately the depths of the potentials in the entrance and exit channels for the same 
reaction, we treat the uncertainties as uncorrelated between the channels involving ^Be and *B, respectively, whereas 
the uncertainties between the two different reactions remain correlated through the use of the same procedure and 
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of the same average values for the renormalization coefficients. When we treated the uncertainties in assessing the 
depths of the potentials in the entrance and exit channels as totally correlated, as we did in Ref . Q , we obtained an 
uncertainty of about 10%. 

We note here that given the observed strong dependence of the calculations on the imaginary part of the potential, 
and the relative independence of the real and imaginary parts of the potential, the measurement of the total reaction 
cross section of ^Be and might be useful. It can set an extra constraint on the renormalization of the imaginary 
part of the potential, and eventually decrease the uncertainty in the potential used, and therefore in the DWBA 
results. 



VII. CONCLUSIONS 

We have measured the elastic scattering of ^Li, ^°B, ^'^C and ^''N on ^Be and ^'^C targets at or around E/A=10 
MeV/nucleon for angular ranges up to around the nuclear rainbow angle, using a fine angle binning of AOiat = 0.5°. 
All these projectile-target combinations, with the exception of -'^'^N+^'^C, have in common the fact that one or both 
partners are weakly bound and we expect contributions from break-up channels to be important. At the same 
time it is known that at energies around 10 MeV/nucleon the contribution of the dynamic polarization potential 
is non-negligible. Parameters for optical model potentials of Woods-Saxon form were obtained from the fit of the 
elastic scattering angular distributions. In addition, nucleus-nucleus potentials were calculated by a double folding 
procedure using six different effective nucleon-nucleon interactions. The nuclear densities calculated for each partner 
in the Hartree-Fock approximation were folded with four different versions of the M3Y nucleon-nucleon interaction 
and with the effective interaction of Jeukenne et al. Q , calculated with two different techniques to account for the 
local density. The resulting nucleus-nucleus potentials were later renormalized to obtain a fit of the elastic scattering 
data. The normalization constants have similar values in all systems for each effective interaction used, which makes 
it appear likely that the procedure can be extended to the calculation of optical potentials for other similar nucleus- 
nucleus systems. From all effective interactions used, we conclude that JLM(l) gives the best results. It provides us 
with an imaginary part that has a geometry which is independent from that of the real part of the potential. The 
imaginary well produced is wider than the real one, as the fit of the data with phenomenological Woods-Saxon wells 
requires. At the same time, it gives the least spread in the value of the renormalization coefficients, which suggests that 
its density dependence accounts very well for the differences between the nuclei involved, particularly in the surface 
region. We find that while the depth of the real potential needs a substantial renormalization ((Ny)=0.366±0.014), 
the imaginary part needs no such renormalization ((Nn/) = 1.000±0.087). This also suggests that the imaginary part 
of the effective interaction is well accounted for. The need for a substantial renormalization of the real part was 
attributed to the effect of the break-up channels, which are very important in nuclei with low binding energies like 
those encountered in our experiments p8[ | . This suggests that the average value of the renormalization constant for 
the real potential depth found above is valid for the region of p-shell nuclei considered in this study and might be 
somewhat different in other regions. 

The renormalized double folded potentials obtained were also used in the DWBA analysis of the proton trans- 
fer reactions with stable nuclei, and the results were found to be in excellent agreement with those given by the 
phenomenological Woods-Saxon potentials. We found that comparison of the double folded potentials and the phe- 
nomenological ones gives a way to select between the different families based on their volume integrals. 

Finally, the procedure found was applied to extract the optical model potentials for the ^Bc and ^B radioactive 
projectiles needed in the description of the ^Be-f ^"B and ^Be-I-'^^'N experiments. Good description for the elastic 
scattering data is found without any need for readjustment of the shape or magnitude of the angular distributions. 
The shape of the angular distributions measured for the proton transfer reactions ^''B(^Be,*B)^Be and ^'*N(^Be,^B)^^C 
are also very well predicted. The calculated DWBA cross sections were used to extract the ANC for the *B^^Be-|-p 
system from each reaction, and consequently the astrophysical 817(0) factor reported in Refs. Furthermore, we 

used the standard deviations of the renormalization coefficients SNy , (5Niy to estimate the contribution of the DWBA 
calculations to the uncertainty of the determined ANC and thus 817(0). We found this contribution to be around 8%. 
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FIG. 1. Spectrum from the elastic scattering of ''Li on the ^Be target at Ei^b = 130 MeV and Oiat = 27.25° ± 0.25°. The 
peak labeled D is a combination of inelastic excitation of ^Be (2.9 MeV) and of double excitation of the target and projectile. 
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FIG. 2. Angular distributions for the elastic scattering of a) (100 MeV) + ^Be, b) "C (130 MeV)+ ^Be, and c) 
(162 MeV)+"C. The curves are fits with the potentials presented in Table II. 

FIG. 3. Same as Figure 2, for the systems a) ^Li(63 MeV)+^Be, b) ^Li(63 MeV)+i^C, c) ■'Li(130 MeV)+^Be, and d) ■'Li(130 
MeV)+^='C. 

FIG. 4. The near-side/far-side decomposition of the elastic scattering for a) ^°B+®Be, b) ^''N+^C, and c) 'Li (63 MeV) + 



FIG. 5. Fit of the angular distributions with the folding potentials of Table IV. The curves are labeled: M3Y/ZR for the 
M3Y zero-range interaction, BDM3Yn for the density-dependent M3Y interactions, and JLM(n) for the interaction of Jeukenne, 
Lejeune and Mahaux, respectively. The cases presented are (a) "N-|-"C, (b) ^°B-|-^Be, (c) ''Li(63 MeV)+"C, and (d) ''Li(63 
MeV)-h^Be. 

FIG. 6. Same as Figure |, but for the systems; (a) ^Li(130 MeV)+"C, (b) "C-f-^Be, (c) ^Li(99 MeV)-f^^C and (d) ''Li(99 
MeV)-h^^Si (data from [||). 

FIG. 7. The renormalization coefficients extracted for the double folding potentials calculated with the six effective nu- 
cleon-nucleon interactions, as described in the text. The projectile-target combinations are those of Table 
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FIG. 8. The double folded potentials calculated with the standard Hartree-Fock mass distributions (dashed lines) are com- 
pared with those obtained when the tail of the proton distribution of ®B is given by the ANC obtained from our experiments (full 
line). Both real (V) and imaginary (W) potentials are shown for the system *B + ^Be, using the JLM(l) effective interaction. 

FIG. 9. The angular distribution for the elastic proton transfer '^*N('^Be,*B)^''C at E;a6=84 MeV, calculated using the 
optical model potential obtained with the JLM(l) effective interaction (dashed line) is compared with the one smoothed by a 
Monte Carlo procedure to account for the experimental resolution (solid line) and with the experimental points. The dotted 
and dash-dotted lines represent the calculated cross section (not smoothed) with the imaginary potential depths renormalized 
by Nw ± SNw (upper panel) . The lower panel presents the ratios of the transfer cross sections calculated using renormalization 
coefficients for the imaginary part of the potential Nvk 4- SNw (dotted line), Nw (solid line) and Nw — Sfiw (dash-dotted line) 
to that calculated with the median value N^. 
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TABLE I. List of the elastic scattering experiments presented in this paper. 



No. 


Projectile-target 


E (MeV) 


^;a5(deg.) 


1 




100 


4 - 30 


2 




130 


4 - 22 


3 




162 


2 - 34 


4 


^Li + ^Be 


63 


4 - 52 


5 


^Li + 


63 


4 - 56 


6 


^Li + 9Be 


130 


4 - 47 


7 


^Li + 13 c 


130 


4 - 47 
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TABLE II. The parameters of the Woods-Saxon optical model potentials extracted from the analysis of the elastic scattering 
data for projectile-target combinations studied here. rc=l fm for all potentials. 



Channel 




Pot. 


V 


W 






■Ay 


aw 








Rv 


Jw 


Rw 








[McV] 


[McV] 


[fm] 


[fm] 


[fm] 


[fm] 






[MeV fm3] 


[fm] 


[McV fm-^] 


[fm] 


(100 MeV) + 


9 Be 


1 


64.2 


30.1 


0.78 


0.99 


0.99 


0.75 


19.8 


1318 


206 


4.51 


136 


4.28 






2 


131.2 


29.7 


0.67 


0.95 


0.90 


0.86 


45.4 


1411 


276 


3.99 


131 


4.46 






3 


203.2 


24.7 


0.81 


1.04 


0.60 


0.83 


61.8 


1428 


499 


3.46 


133 


4.59 


"N (162 MeV) + 


13c 


1 


79.22 


30.27 


0.96 


1.05 


0.76 


0.72 


17.4 


1542 


221 


4.52 


105 


4.69 






2 


134.76 


35.23 


0.88 


1.05 


0.75 


0.67 


18.3 


1525 


299 


4.28 


120 


4.61 






3 


176.03 


35.84 


0.86 


1.07 


0.72 


0.65 


23.3 


1527 


361 


4.15 


125 


4.62 






4 


241.36 


37.45 


0.82 


1.06 


0.71 


0.66 


27.5 


1533 


438 


4.00 


129 


4.61 






5 


306.44 


39.14 


0.81 


1.05 


0.68 


0.68 


36.1 


1552 


522 


3.90 


132 


4.61 


13c (130 MeV) + 


«Bc 


1 


94.2 


20.9 


0.77 


0.99 


0.87 


0.97 


15.0 


1592 


223 


4.19 


94 


4.96 






2 


164.2 


23.0 


0.67 


0.98 


0.86 


0.95 


14.2 


1576 


283 


3.94 


99 


4.87 






3 


226.7 


31.9 


0.62 


0.90 


0.85 


0.95 


14.8 


1573 


328 


3.81 


113 


4.70 


■^Li (63 MeV) + 


9Bc 


1 


134.4 


19.82 


0.54 


1.03 


0.95 


0.92 


8.0 


1414 


267 


3.90 


137 


4.66 






2 


221.6 


27.33 


0.54 


0.92 


0.83 


0.97 


10.6 


1449 


367 


3.50 


153 


4.60 






3 


276.5 


37.3 


0.61 


0.81 


0.72 


1.02 


15.7 


1482 


499 


3.27 


158 


4.64 


'Li (63 MeV) + 


Uc 


1 


54.3 


29.9 


0.92 


1.03 


0.79 


0.69 


28.8 


1318 


209 


4.21 


144 


4.26 






2 


99.8 


22.0 


0.77 


1.01 


0.81 


0.81 


21.6 


1363 


257 


3.92 


109 


4.49 






3 


154.8 


22.7 


0.76 


1.00 


0.71 


0.83 


19.8 


1378 


357 


3.64 


111 


4.51 






4 


244.6 


26.4 


0.68 


0.96 


0.71 


0.84 


20.4 


1382 


438 


3.47 


117 


4.45 


'Li (130 MeV) + 


^Be 


1 


60.0 


17.71 


0.86 


1.07 


0.65 


1.12 


150. 


1564 


217 


3.58 


154 


5.33 






2 


129.4 


30.7 


0.57 


0.80 


0.90 


1.17 


208 


1488 


261 


3.77 


158 


5.02 


'Li (130 MeV) + 




1 


123.3 


32.74 


0.76 


0.94 


0.76 


0.90 


79.1 


1406 


297 


3.79 


145 


4.66 






2 


157.9 


31.97 


0.63 


0.90 


0.87 


0.94 


77.3 


1393 


289 


3.83 


133 


4.59 






3 


201.9 


25.59 


0.73 


1.03 


0.69 


0.86 


129. 


1418 


419 


3.52 


142 


4.66 






4 


300.0 


30.78 


0.73 


0.98 


0.64 


0.89 


147. 


1441 


543 


3.37 


150 


4.63 
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TABLE III. Radii and binding energies of the calculated Hartree-Fock one-body densities, compared with the experimental 
data. Rp, Rn, Rm and Hch stand for the root mean square radii of the calculated proton, neutron, mass and charge distributions 
respectively and R^^'' is the experimental charge rms. B are the binding energies. 



nucleus 


R 

[fm] 


R 

[fm] 


R 

[fm] 


R 

[fm'i 


[f 


[MeV] 


[MeYL 
Ref. jij 


Rcf. fci 


6Li 


2.21 


2.20 


2.20 


2.21 


2.54(5) 


31.929 


31.994 


^Li 


2.15 


2.35 


2.26 


2.16 


2.39(3) 


39.234 


39.244 


■^Be 


2.37 


2.14 


2.28 


2.38 


2.36(2)* 


37.606 


37.600 


8B 


2.57 


2.18 


2.43 


2.58 


2.45(5)* 


37.744 


37.737 


9Be 


2.26 


2.39 


2.33 


2.29 


2.50(9) 


58.203 


58.164 


10b 


2.40 


2.39 


2.40 


2.45 


2.45(12) 


64.631 


64.750 


12 c 


2.44 


2.43 


2.44 


2.49 


2.47(2) 


92.149 


92.161 


13 Q 


2.47 


2.56 


2.52 


2.53 


2.440(25) 


97.135 


97.108 


14n 


2.58 


2.57 


2.57 


2.64 


2.58(2) 


104.246 


104.658 



* proton density rms radius obtained by Tanihata at al P5| from interaction cross sections. 
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TABLE IV. Best fit renormalisation paxameters Nv and Nw for folding potentials with various effective interactions (see 
eqs. (15) and (16)). For each reaction channel, the values of Nv are given in the first line and Nw in the second line. For 
each effective interaction, the mean values and dispersions are given in the last two lines. Only cases 2-7 are used to determine 

averages, as described in the text. 



No. 


Projectile-target 


JLM(l) 


JLM(2) 


M3Y/ZR 


M3Y/FR 


BDM3Y1 


BDM3Y3 


1 


i^'N (162 McV) + ^^C 


0.4-56 


0.509 


0.778 


1.275 


0.721 


0.832 






0.844 


0.996 


0.469 


0.887 


0.419 


0.492 


2 


^"B (100 MeV) + ^Be 


0.368 


0.387 


0.516 


0.667 


0.584 


0.668 






1.168 


1.131 


0.571 


1.116 


0.506 


0.596 


3 


^■^C (130 MeV) -1- ^Be 


0.369 


0.413 


0.489 




0.576 


0.648 






0.937 


1.124 


0.726 




0.550 


0.656 


4 


^Li (63 MeV) + ^''C 


0.323 


0.364 


0.588 


0.787 


0.552 


0.634 






1.00 


1.007 


0.503 


0.831 


0.458 


0.535 


5 


'Li (63 MeV) + ^Be 


0.360 


0.403 


0.588 


0.759 


0.568 


0.645 






1.00 


1.438 


0.818 


1.175 


0.733 


0.864 


6 


'Li (130 MeV) + '-^C 


0.380 


0.418 


0.595 


0.914 


0.571 


0.651 






0.957 


1.077 


0.508 


0.893 


0.472 


0.547 


7 


Li (130 MeV) + Be 


0.368 


0.413 


0.489 


0.806 


0.576 


0.648 






0.937 


1.124 


0.726 


1.110 


0.550 


0.656 


8 


"Li (99 MeV) + ''^C 


0.449 


0.493 


0.716 


1.178 


0.687 


0.785 






1.044 


1.166 


0.536 


0.942 


0.510 


0.585 


9 


"Li (99 MeV) + ^»Si 


0.368 


0.408 


0.565 


0.960 


0.534 


0.611 






1.168 


1.324 


0.683 


1.170 


0.621 


0.726 


10 


'Li (63 MeV) + ^^C 


0.278 


0.309 


0.502 




0.478 


0.546 






0.746 


0.920 


0.464 




0.423 


0.493 


11 


'Li (79 MeV) + ''^C 


0.315 


0.347 


0.521 




0.505 


0.573 






0.864 


1.009 


0.458 




0.426 


0.493 




average of cases 2-7 


0.366±0.014 


0.405±0.017 


0.553±0.062 


0.787±0.089 


0.578±0.010 


0.658±0.013 






1.000±0.087 


1.143±0.145 


0.631±0.131 


1.025±0.153 


0.553±0.082 


0.631±0.115 



TABLE V. Parameters of volume Woods-Saxon type potentials that best fit the nuclear part of the numerical poten- 
tials obtained with the double folding procedure using the JLM(l) effective interaction in the range r=2-12 fm (see text). 
Renormalization of the depths is included. Kv and Hw are the half-radii of the potentials. 



Projectile-target 




V 


W 




Rw 


ay 




Jy 






[MeV] 


[MeV] 


[MeV] 


[fm] 


[fm] 


[fm] 


[fm] 


[MeV fm3] 


[MeV fm3] 


'Be + ^"B 


84 


63.8 


29.4 


3.18 


3.49 


0.85 


0.95 


210 


130 


^B + ^Bc 


81 


67.0 


31.8 


3.18 


3.54 


0.88 


0.99 


236 


145 


^Bc + 


84 


79.1 


36.0 


3.30 


3.62 


0.88 


0.98 


207 


126 




78 


85.2 


39.3 


3.30 


3.76 


0.91 


1.02 


216 


145 
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